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Abstract. We study the dynamics of soliton solutions to the perturbed mKdV 
equation dtu = d.^{~dlu - 2u^) + eVu, where V G Cl{M), < e < 1. This type 
of perturbation is non-Hamihonian. Nevertheless, via syniplectic considerations, 
we show that solutions remain 0(e(t)^/^) close to a soliton on an 0(e^^) time 
scale. Furthermore, we show that the soliton parameters can be chosen to evolve 
according to specific exact ODEs on the shorter, but still dynamically relevant, time 
scale 0{e~^^'^). Over this time scale, the perturbation can impart an 0(1) influence 
on the soliton position. 



1. Introdution 

We consider the modified Korteweg-de Vries (mKdV) equation with a small external 
potential 

(1.1) dtu = d^{-dlu - 2m^) + eVu . 

where < e ^ 1, G C^(M), i.e. V and V are continuous and bounded. 



The unperturbed case of (1.1), 



(1.2) dtu = d^,{-dlu-2u') 

is globally well-posed in for /c > 1 (see Kenig-Ponce-Vega [H]), and possesses single 
soliton solutions u{x,t) = rjlxja + c^tjc), for a G M and c G M \ {0}, where r]{x,a,c) = 
cQ{c{x — a)) with Q{x) = sech(x) (so that —Q + Q" + 2Q^ = 0). The solitons are 



orbitally stable as solutions to the unperturbed mKdV (1.2) (see [31 HI [28], [7]), i.e. 
the solutions stay close to the soliton manifold 



M = { r]{x, a,c) \ a eR,c> 0} 



if they are initially close. 



Our first main result, Theorem 1.1 , shows that this type of orbital stability remains 



true for the structurally perturbed mKdV ( 1.1 ), in the following sense: solutions which 



start an distance u from the soliton manifold M remain within an distance 



1.2 



{tu + et^/^)e'"^* up to time e~^loge"^. Our second main result result. Theorem 
shows that on the shorter time scale e~^/^loge"^, we can predict the location on 
the soliton manifold by solving a system of two ODE for the position parameter a 
and scale parameter c. Strong agreement between this prediction and the numerical 
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solution of (1.1) is illustrated in Fig. |1.1| and Fig. 1.2 We prove the global well 
posedness of (1.1) in H^, by adapting the argument of Kenig-Ponce-Vega 
Apx. |X| 

The forced KdV equation 



(1.3) 



dtu = d^i 



3u') + ef 



is a model for free-surface shallow water flow [20] with contributions to / arising from 
surface pressure and bottom topography. Numerics and experiments discussed in pU] 
show that this type of perturbation can effect the evolution of a single soliton by 
generating a procession of small solitons ahead of, and dispersive waves behind, the 
primary soliton. 



Both (1.1) and (1.3) are specific instances of a family of gKdV equations with 
general perturbation 

dtU = d^{-u^^ - u'P) + ef 

for ^ G N, p > 2, and / = f{x,t,u). The case p = 3 (mKdV) is the unique member 
of the gKdV family that avoids a certain anomaly with the symplectic structure. 
Specifically, for p = 3, one has d'^dcTj G but this fails for p ^ 3. For p = 3, one 
can symplectically project onto the tangent space of the soliton manifold M rather 
than on a skew space. The difference between p = 3 and p 7^ 3 is illustrated in the fact 
that the local virial estimate of Martel-Merle ||21j simplifies for p = 3. Nevertheless, 
we believe that the analysis of the paper carries over in some form to p 7^ 3 and more 
general / of the form f{x,t,u). We chose (1.1) as the mathematically simplest case 
in which to illustrate our method. 

1.1. Statements of main results. 

Theorem 1.1 (orbital stability). Let 6 > and oq, cq G M such that 25 < cq < {26)~^ . 



Suppose u{x,t) solves (1.1) with initial data u{x,0) such that 



u = ||m(x,0) -r7(x,ao,Co)||//i < e^/^ 

Then there exist trajectories a{t) and c{t) so that the following hold, where T is 
the maximum time such that 6 < c{t) < 6^^ for all < t < T and w{x,t) == 
u{x,t) — ri{x,a{t),c(t)). First, we have the following bounds on the deviation w: 



(1.4) 



Second, we have T > C ^ and the following estimates for the trajectories a{t) and 
c{t): 
(1.5) 



ec ^{Vt], (x 



The constants C in (1.4), (1.5) depend on ||V^||ci CL^d 5. 
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We remark that the same result holds for cq < 0, since ri{x, a, —c) = —ri{x, a, c). 



(1.6) 



Theorem 1.2 (exact predictive dynamics). Suppose u{x,t) solves (1.1) with initial 
data u{x, 0) satisfying 

u = \\u{x,0) - r]{x,ao,Co)\\Hi < e^/^ 
where Cq > 0. Let {a{t),c{t)) evolve according to the ODE system 

a = (? + ec~^{Vrj, {x — a)ri) 
c = e{Vrj, rf) 

with initial data a(0) = a^, c(0) = Cq. Then for 

0<t<T = ae^^/^loge"^ , a = a(co, ||l^||ci) > 0, 
we have the following estimates with w{x, t) = u{x, t) — ?7(x, a(t), c(t)) : 
(1.7) \\w\\l^ m + lle-^l'-'^lMiLa < C(uj + et^/^)e^'"^'' . 

where C = C{cq, \\V\\ci). 
We remark that if one selects initial data so that u < e^/'^, then the two terms on 



the right-side of the estimate (1.7) balance on the e time scale. In this case the 
bound becomes e^/^e'"'^^''^*. 

1.2. Relation to recent work. The energy-Lyapunov based methods for proving 
orbital stability of solitons subject to perturbations (of the data, as opposed to the 
structural perturbations considered here) were developed by Benjamin [3], Bona [1], 
Weinstein [28], Grillakis-Shatah-Strauss [HI |T2]. In the last decade several results 
have emerged using the same basic framework to address the dynamics of solitons 
for equations subject to structural perturbations PIHllIOlIISllHIISlIIZlIIHllHlIIlEl 
[231 El]. The nonlinear Schrodinger equation (NLS) with slowly varying potential was 
considered by Frohlich-Gustafson-Jonsson-Sigal [9] and a result of "orbital stability" 
type was obtained, however the estimates were not strong enough to obtain "exact 
predictive dynamics" . Holmer-Zworski [18] obtained exact predictive dynamics plus 
refined accuracy by adopting the conceptual perspective of symplectic projection, but 
also, at the technical level, finding an appropriate distortion of the soliton manifold 
that enabled refined Lyapunov estimates. This "symplectic projection plus correction 
term method" has been subsequently pursued in different contexts in Datchev- Ventura 
[S], Holmer-Lin [13], Holmer-Perelman-Zworski [IB], and Pocovnicu [2S]- To treat a 
problem in which the perturbation gives rise to significant dispersive radiation, a 
different approach was employed by Holmer [13]. He treated the KdV equation with 
a slowly varying potential, and used the Martel-Merle local virial estimate [23] |23] 
to supplement the energy Lyapunov estimate. In this paper, we follow this approach 
as well. We show the method is sufficiently robust to handle small non-Hamiltonian 
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Figure 1.1. With external potential given by Vi as in (1.8), the top 
plot gives the rescaled evolution U{X, T), the bottom two plots give the 
comparison between the evolution of the parameters obtained solving 
the ODE system and exact PDE evolution, i.e. we fit the solution to 
ri{X, A, C), and plot T versus A and C respectively. 



perturbations, which had not been considered in any of the above papers. A stochastic 
variant of the problem we consider has been addressed by de Bouard-Debussche [5] 
without the use of the local virial estimate. Work in progress by Holmer-Setayeshgar 
[T^ will adapt the present paper to the stochastic setting and obtain a refinement of 
the results of [5]. 



1.3. Numerics. To solve (1.1) numerically we adapt the method in [26] which is 
based on the fast fourier transform in x, then fourth-order Runge-Kutta for the re- 
sulting ODE in t. We use the rescaled coordinate frame X = e~^^^x, T = e~^t, and 
consider the equation on [— 7r,7r). If U{X,T) solves 

OtU = dx{-dlU - 2f/3) + V{X)U , 
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Figure 1.2. These plots are analogs of Fig. the external potential 



is given by V2 as in (1.9). 



with initial data 

f/(0, X) = r]{X, Ao, Co) = e^/^ao, e-'/%) , 
thenM(a;,t) = e^^'^U{e^^^x, et) gives a solution of (1.1) on [— vr/e"^/'^, vr/e^/^) with initial 



data u{0,x) = ■/^(x, ao, cq), and periodic boundary conditions. Fig. |1.1| and Fig 1.2 
plot the evolution of the soliton initial data (after rescaling) in the following external 
potential respectively 



(1.8) 
(1.9) 



Vi = -10 cos2(6X) + 6 sin(lOX) 
V2 = 8cos^(4X) -4sin(2X). 



Note that to examine the solution u{x,t) on time interval < t < Ce ^/^(or Ce ^) 
we should let U{X,T) evolve for time Ce^''^(or C). 
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2. Background on Hamiltonian structure 

Let J = dx, and consider L^(]R i— ?► M) as a manifold with metric (fi, V2) = J V1V2 dx, 
we can define the symplectic form as 

(2.1) Uj{Vi, V2) = {ViJ'^V2) = {Vi, d^^V2) , 

where is given by 



dcf 



+00 



f{.y) dy. 



The mKdV equation (1.2) is the Hamiltonian flow associated with 

1 
2 



Hoiu) 



i.e. we can write (1.2) as 
(2.2) 



dtu = JH'oiu) 



Solutions to mKdV also satisfy conservation of mass M{u) and momentum P{u) 
where 

M{u) = I udx, P{u) = ^ 



u dx . 



We define 2-dimensional manifold of solitons M as 

M = { r]{-, a,c) \ a eR,ceR\ {0}} . 



The symplectic form (2.1) restricted to M is given by u\m = da A dc. We denote 
T] = ri{-,a,c), the dependence of (a, c) on t] is always meant implicitly. The tangent 
space at rj is given by 

Tr/M = span{ daf], d^r] } . 



Note that JH^ir]) G T^^M, thus the flow associated to (1.2) will remain on M if it is 
initially. Specifically, direct computation shows 



(2.3) 



JH',{7])=c''daV. 



which, together with (2.2), explains the form of the expression for single solitons. 



This is equivalent to saying that the flow (2.2) restricted to M (and thus stays on 
M) is given by 

d = 

(2-4) , . „ 

c = 
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One can also get (2.4) by first restricting Hq to M to obtain 



and then noticing that (|2.4|) is just the solution to the Hamilton equations of motion 

\m- 



for Hq{vi) with respect to uj\ 
(2.5) 



a = — -— = c 
oc 

c=^ = 
da 



Note that we can write (2.3) as 

(2.6) JH'o{r]) + c^JP\r]) = 0. 
From this, we learned that L^t]) = 0, where 

(2.7) L{u) = Ho{u) + c^P{u) . 

which is the Lyapunov functional used in the classical orbital stability theory, see 
[28]. 

Next, we define the symplectic orthogonal projection operator at (a,c): 
by requiring that 

(ni J, j-'dav) = (ni J, j-'d,v) = , 

where U-^,, = I — Ila,c, equivalently, 

Ha.c/ = (/, J~'^dcr])daV - if, J~^dar])dcr] . 

Note that for mKdV, 

J^^daTj = —1] and J'^dj] = c~^{x — a)rj . 

3. Decomposition of the flow 
We can arrange the modulation parameters a{t) and c{t) so that 

^a{t),c{t) [u{x, t) - ?7(x, a(t), c(t))] = . 
This is a standard fact and we recall it in the following 

Lemma 3.1. Given a, c, there exist 5i > 0, C > 0, such that if u = ri{-,d,c) + w 
with \\w\\fji < 6i, then there exist unique a, c such that 

(3.1) w{x, t) =^ u{x, t) — ri{x, a{t), c{t)) 
satisfies the symplectic orthogonality conditions 

(3.2) {w, J-^daV) = {w, J-^dcV) = . 
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Moreover, 

\a-a\<C\\w\\Hi, |c-c| < . 
Proof. Define : if^ x M x M ^ 

(p{v,a,c) = 

\_{v ~ ri,[x — a)ri) 

Using oj\m = da A dc, we can get the Jacobian matrix of (f) with respect to (a, c) at 
('?(■, a, c), a, c) 

1" 



(^a,c0)(^?(-,a,c),a,c) 



1 



which imphes, by the imphcit function theorem, that the equation a, c) = can 
be solved for (a, c) in terms of m in a neighborhood of ri{-, d,c). □ 



Now since u = w + rj and u solves (1.1), we compute 

dtw = dx{—d1w — Qrfw — Grjw^ — 2w^) + eVw — Fq 



(3.3) 
where 



dx{Cw — c^w — Qriw"^ — 2w^) + eVw — Fq . 



and Fq results from the perturbation and dt landing on the parameters: 

Fo = (d - c^)dari + cdcV - eVr] . 

Next, decompose Fq into the symplectically parallel part Ilac-Po and symplectically 
orthogonal part II^^Fo, explicitly, 

(3.4) n,,,Fo = {d-c'- e{Vv, J-'d,v))daV + (c + e{Vv, J-'daV))dcV , 

(3.5) ni.Fo = -eVr] + e{Vr], r'd,r])daV - e{V7], J~'daV)dcV • 
We now obtain the equations for the parameters: 



Lemma 3.2 (effective dynamics). Given V G C^, suppose that w defined by (3.1) 



satisfies the orthogonality conditions (|3.2|). Then there exists a > such that 
(3.6) 



\dtr] - c^daV - eUa,ciVr])\\T, 



cM ^ 



+ e e 



-a|x— a| 



Explicitly, 



(3.7) 



\d-c'-e{Vr],J-'dcV)\ < ||e-"l^-'^lu;||^i+e||e-'^l^-'^l«;||Hi, 
|c + e{Vr], J~'daV)\ < We-'^^^'Mm + e||e-"l^-'^lw||^,i . 



As all norms on a finite dimensional space are equivalent, we can take 
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Proof. Recall that 

dtw = JH'^{ri)w - J{Qriw'^ - 2w^) + eVw - Fq . 



Write TZ for the error terms of the same order as the right hand side of (3.6), take 
derivative with respect to t for {w, J'^dat]), we have 

= {dtW, J-^daV) + {W, J-'dadtV) 

= - {Fo, J~'dav) + {JH'^{v)w, J-'dav) + {w, J-'dadtv) + 7^ 
= - {Fo, J-'dav) + {w, J~'daidtv - JKiv))) + ^ 

= - (Fo, j-'daT]) + {w, j-i9,(n,,,Fo)) + 7^ , 

where for the penultimate equality we have used J* J^^ = —I and the self-adjointness 
of Hq, and for the last that 

dtv - JKir]) = {a- c')daV + cd,v = n,,,Fo + 0{e)daV + 0{e)d,v ■ 

Taking derivative for {w, J^^dct]), similar computation gives 

= -{Fo,r'd,r]) + {w,J-'d,{Ua,cFo))+n. 



Combining with (3.8), and applying the orthogonality conditions for the second terms 
when da and dc land on the coefficients of Ua^cFo, the lemma follows from Cauchy- 
Schwarz and the smallness of w. □ 

4. Local virial estimate 

In this section we review, and then apply, part of the local virial estimates due to 
Martel-Merle. Let $ G C{R), ^(x) = $(-a;), $' < on (0, oo), such that 

<^{x) = 1 on [0, 1] , ^{x) = e"^ on [2, oo) , e"^' < $(x) < 3e"^' on [0, oo) . 

Let '^{x) = $(?/) dy, and for A 0, set "^Ai^c) = A'^{x/A), we have following 

Lemma 4.1 (Martel-Merle f2T\ [22] local virial spectral estimate). There exists A 
sufficiently large and Aq sufficiently small, such that if w satisfies the orthogonal 



condition (3.2), then 

-{'^Aix-a)w,d^Cw) >\oj {wl + w'^)e-\''-^\'^dx. 

Denoting ■?/;(■) for ^ a{' ~ c^), we now proceed as in [21]: 
Lemma 4.2 (local virial estimate). Suppose V is hounded, then there exist a > and 



Kj > 0, j = 1,2, such that if w solves (3.3) and satisfies the orthogonality conditions 



(3.2), then 



(4.1) lle""'''""'^!!^! < -Kidt J 'ipw^ dx + K2e^ + K2e||w||^i . 
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Proof. From the equation for dtW, we have 
dt [ ipw"^ = — a I ip'w^ + 2 /" ipwdtw 

I.I .."2 



= - d y i)'w^ + 2 J ^wd^{Cw) 

— 2(? J ipwdxW — 12 J ipwdx{i]w'^) 

-4 J ijwdxiw^) + 2e j ijVw^ - 2 / tpiuFo 
Using integration by parts, 



^ I + II 



^ III + IV 



^ V + VI + VII 



hence 
(4.2) 



11 + 1111 



[a — c 



111 = 



1p W 



4w\\l. + \\e-''^''-Mm\H\Hi 



by (|3.7|). Following from the boundedness of ip and V, and the estimate 

|IV| < ||e-"l^-"li 



< 



(4.3) 



we obtain 

Itt7-I ^ II —rfln^ — n, ,, ^ ,, ,, 

|V| = |3 y"^V|<||u;||^.||e-l^-^l/(2-^)^lli., 
|VI|<e||«;||^,, 

where for the second estimate we have used ip' = $((x — a) /A) and the definition of 
$. Decomposing VII term as 

VII = -2 1 t/jwUFo -2 1 ipwU^Fo = VIIA + VIIB , 



we have by Lemma 3.2 that 



(4.4) 



VIIA < e\\w\ 



HI 



\e-''\^-Mh\H\m, 



and by II-'-Fo ~ et] (see (3.5)) that for any /i > 0, 

(4.5) VIIB < e||e-"l^-'^l^||Hi < /i^'e^ + /i||e'"l^-'^l^||^i 

Note in above estimates the value of a may change from one line to the next, but we 
can choose one single small enough a that works for all. 



By Lemma 4.1, we have 



II = 2{i;w, dx{Cw)) < -Ao j {wl + w2)e-l"-"l/^ dx . 



Combining with (4.2), (4.3), ( |4.4 ) and (4.5), the estimate (4.1) follows by the smallness 
of taking A large enough so that 1/(2^4) < a, and > suitably small. □ 
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In this section we formulate the energy estimate necessary for the estimation of the 
error term w. Recall C = —d^ — Qrf + c^. Let 

£ = - {Cw, w) — 2 I rjw^ dx — - f w'^ dx , 



Note that C = H'^^iri) + = L"{ri), see (2.6) and (2.7). We have classical coercivity 
properties for £ (for a proof, see e.g. Prop 2.9] or [T71 Prop 4.1] for a more direct 
proof - note that £ is the operator L+ considered there): 

Lemma 5.1 (energy spectral estimate). Suppose that w satisfies the orthogonality 



condition (3.2). Then 
(5.1) 



{Cw,w) > \\w^ 



x\\l^ 



+ c'\w\ 



L2 , 



Since we impose a lower bound on c in Theorems |1.1| it follows from (5.1) that if 
is smaller than some (e independent) constant, then 

Mt)\\l.-£{t) 

Lemma 5.2 (energy estimate). Suppose we are given V G Cl, Sq > and w{x,t), 
such that 5o < c{t) < 5q^ , w solves (3.3) and satisfies the orthogonality conditions 
(g, then 

(5.2) m < e\\w\\l. + e||e-"l— + Ml.Je-'^^^-Mm + \H\%i ■ 

where the implicit constant depends on 6q, (Tq and the bounds on V and V . 

Proof. We compute 

dt£ ={Cw, dfw) + cc||w||^2 — 6{{ddar] + cdcr])r]w, w) <(— I + II + III 

- {dtw, 6r]w^ + 2w^) - 2{{ddaV + cdcV),w^) ^ IV + V 

Substitute ( 3.3[ ) into I: 

I = {Cw, d^{Cw)) - c^Cw, d^w) - {Cw, d^{6r]w'^ + 2w^)) + {Cw, eVw) - {Cw, Fq) 
= lA + IB + IC + ID + IE 

First, lA = 0. Integration by parts yields IB = ~6c'^{r]r]x,w'^). By the boundedness 
of V and V, 

IB < e\\w\\l,^ , 

and since C{TM) C TM (by direct computation), we have 

IE = -{Cw, UFo) - {Cw, n^Fo) = -{Cw, U^Fq) , 



but by (3.5) 



\{Cw,n^Fo)\<4e-''^''-Mm 



12 
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|IE| < e||e-"l"-'^lM;||j^i . 

Combining, we obtain 

(5.3) I = IB + IC + O (elkll^i + e||e-"l^-"lw;||jyi) 

= -Qc^i-qr]^, w'^) - {Cw, d^{6riw^ + 2w^)) + O (e||w;||^i + e||e-"l^-°lw;||H3 
Substituting ( 3.3[ ) into IV, we have 

(5.4) IC + IV = - {d^i-c^w - 6r]w^ - 2w;^) + eVw - Fq , 6r]w^ + 2w^) . 



By (3.7), we have 



(5.5) 
hence 

Note 



\a-c'\<e+\\e-"^'^-Mm, 



\{Fo,6r]w^ + 2w^)\<e\\w\\l^ 



icl <e+||e-"l"-"lw||ii 



Hi 5 



\W\ 



-ok-a|^||2 



HI 



■{d^{-c^w), Qrjw^ + 2w^) = -2(? / r^V dx . 



Estimating the rest of the terms in (5.4) using Cauchy-Schwarz and that 
we obtain 

(5.6) IC + IV= -2c2(r/',w^) +0(e||w||^i + He" 



< 



By (5.5) again, and that drji = —dai], we have 
(5.7) 11 + V = 2d(r/',w3) + 0(e||w||^i + ||e 

and 

IB + III =6(a - c^)(vVx, w^) - Q(c(dcV)v, w"^) 

(5.8) 



e\\w\\l. + \\e-''\''-M\m\\w\\l.. 



Apply (5.5) again to the sum of (5.6) and (5.7), then combine with (5.3) and (5.8) 



we can obtain (5.2). 



□ 



6. Proof of the main theorems 



First, we give the proof of Theorem 1.1 



Let [0, T'\ be the maximal time interval so that 
(6.1) 



for /i > chosen small enough to ensure the validity of Lemmas 3.1, 3.2 |4.2[ and 



5.2, and also small enough to beat some constants in the estimates that follow (as 
explained below). 
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Let 



dcf 



-a\x~a{s)\ 



wis)\\jji ds 



J^{t) = sup 

0<s<t 



Integrating the local virial estimate (4.1) gives 
(6.2) 



< J^it) + eH + e / T{s) ds . 



Integrating (5.2) over < t < r yields 



S{r) < S{0) + e / J^{s) ds + eri/2v(r) + -F(r)V(r) + rj^(r)^ . 



Using that S{t) ~ II'^(t) and then taking the sup of the above estimate over 
< r < t, we obtain 



Jo 

By (Q and the estimate et^/^V{ty/^ < fi-^H + ^iVit) this implies 

^(t) ^ e f J^(s) ds + J^(0) + ^'^H + /iV(t) 
io 

Substituting (6.2) into here, taking /i (introduced in (6.1) above) small enough to 
beat the implicit constants, 



(6.3) 

Hence, for some k > 0, 

d 
dt 

Integrating yields 



J^{t) < e / J^{s) ds + T{0) + eH . 



T{s)dsj < e-^'\T{0) + eH) 

fj^{s)ds<{e^^'-l){e-'j^{0) + l) 
Jo 



Substituting this back into (6.3), 

J^{t) < e^'^'TiO) + e((e"^* - 1) + et) 
For the second term, we might as well bound (e'*'^* — 1) + et < ete"^*, so 

This enables us to reach time ae~^ loge~^, for a > small, while still reinforcing the 
bootstrap assumption (6.1). Returning to (6.2), we obtain the bound for V(t), thus 
completing the proof of (1.4). The -^^[ot] estimates (1.5) follow from integrating (3.7) 
in time and applying (1.4). The L'^j^-^ estimates also follow from (3.7) by dropping 
the spatial localization in the terms on the right-hand side of (3.7) and applying the 
bound on m given by(1.4). 

[0,T] ^ 1 



14 



QUANHUI LIN 



Now we discuss the proof of Theorem 1.2 
Let a, c solve the ODE system 

a — — ec^"^(K?7, (x — a)r/) = 

c — €{Vfi, fj) = 

with initial data a(0) = ao, c(0) = cq, where = cQ{c{x — a)). Since |c|, |c| < e, we 
can assume Sq < c, c < 6q^ on [0,T]. Define 

a = a — a , c = c — c, 

we have 

{Vt], {x - a)r]) - {Vfj, {x - a)f/) = j3i{a - a) + /32(c - c) , 
where we have defined 

/3i = f (v(- + a)-V(- + d)) xif dx , 

a — a J \ c c / 

(32 = —. [ (vi- + a)-V{^ + a))xr]'dx, 
c — c J \ c c / 



similarly, 



where 



-{Vri.ri) - l{Vri,ri) = 7i(a - a) +72(0- c) 
c c 



71 



72 



/" fv(? + a)-r(? + a))r/2dx, 
. — a J \ c c / 

/" fv(- + a)-y(? + a))r/2dx, 
— c J \ c c / 



Denote T^i, 7^2 for the error terms in Lemma [3.2[ i.e. 

d — — ec^^{Vri, (x — a)ri) — T^i = 
c-e(Vr7,r7) -7^2 = 0, 
Apply (1.5) to (3.7), we obtain 
(6.4) \\njhi^,<Citu + ee/ye'''"'"' 

Note 



J = l,2. 



and 

denoting 
and 



c c c c _ 
c c c cc ' 

cd — cd = cd + (c — c)d . 



61 = - [{c^ + c^ + cc) - (c' + ec~\Vfi, {x - d)f,)) + e/32] 
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(6.5) 
Writing 



a 


1 


c 





ec7i €{02 + C72) 



+ 



7^l 
7^2 



Ait) 



ec7i 6(6*2 + C72) 



From the boundedness of /3j, 7^, 6j, j = 1, 2, which is a resuh of the boundedness of 
V, V, c and c, we have the estimate 



(6.6) 



< 



e 1 
e e 



Writing p{s) = [ea? + c^)^^^, then by above estimate 

\P\ ^ - WMa\ + |c| + |7^l|) + |c|(e|a| + e|c| + |7^2|)] 

< - re(ea2 + c^) + e^/^(ea2 + c^) + e|a||7^l| + |c||7^2|l 

<el/^p + eV2|7^l| + |7^2|. 
By Gronwall and p(0) = 0, we obtain 

p{t) < Ce""'"'"' f (el/2|7^l| + |7^2|) (s) ds . 
Jo 



Applying (6.4), we obtain 



pit) < Ce^'"'^\uj + tt 



AI2V1 



recalhng the bounds on t and uj in Theorem 1.2, this gives 

pit)<Ct^l\uj^et^l^)e''^'^"\ 

The bounds on a and c now follow from the definition of p: 

|a| <C(a; + eti/2)gCeV2,^ 



Compare the above two estimates with (1.7), we can conclude the proof of Theorem 

Remark 6.1. The e~^/^ constraint on the time scale stems from the fact that the 

are only of order e^^^. 



eigenvalues of 



e e 
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Appendix A. Local and global well-posedness 

The global well-posedness for gKdV in energy space was obtained by Kenig- Ponce- 
Vega in [in], where they introduced new and powerful local smoothing and maximal 



function estimates, especially, they proved the local well-posedness for (1.2) in H 



for s > 1/4. To prove well-posedness for (1.1) at level of regularity, the full 
strength of these estimates is not needed, we here follow the presentation of [16] Apx. 
A and make necessary modifications. 

Let Qn = [n — + and Qn = [n — 1, n + 1]. An example of notation is: 



Note that due to the finite incidence of overlap, we have 



\U\\fooT2 t2 ~ \\U\\fooT2T2 



We omit the e in (1.1), and consider 



(A.l) dtu = d^{-dlu - 2u^) + Vu, V eCl. 



As in [16] , we first prove a local smoothing estimate and a maximal function estimate 
(weak versions), by an integrating factor method: 

Lemma A.l. Suppose that 

(A. 2) vt + v^^^ -Vv = f , 

then there exists C > 0, such that if 

T<C{1+\\V\\l^)-\ 
we have the energy and local smoothing estimates 



and the maximal function estimate 

The implicit constants are independent of V . 



Proof. Let (j){x) = — tan ^{x — n), and set w{x,t) = e'^^^^v{x,t). By (A. 2), 

dtw + Wxxx - 30'w7,, + 3(-0" + i(l)')^)w, + (-0'" + 30>' - i(t>'f)w -Vw = f 
integrating its product with over x, 

dt\\w\\l2 = -6(0', wl) + {-r + 2(0')^ nj') + m + 2(e'^/, w) , 
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integrating this identity over [0,T], and using 
\\wiT)\\l2+ 6\\{x - n)-^w^\\l2r2 



{x — n) , we obtain 



dt . 









J e'^fw dx 


Jo 





for some constant Ci > 0, replace T by t, and take supremum over t G [0,T], we 



obtain, for T < ^{l + \\v\\l°°), the estimate 



\W\\ TOOT2 



+ ll(x-n) ^w4l2r2 < IkoIlL + 









J e'^fw dx 


^0 





dt . 



note that < e"^/^ < ^^x) < ^ 

oo, we can convert the above estimate back to 

an estimate for v: 

-T 



I^0||l2 + 



^^'^fv dx 



dt. 



Estimating as 



e^'t'fv dx 



dt < 



l}^li\\v\\l^li , 



and then taking the supremum in n yields the second estimate in (A.3). Estimating 
instead as 

{d-^fd^{e^^v)dx dt 











f 


j e^'^fvdx 


dt= 1 


1 






Jo 





t2 t2 



poo r2 t2 



and then taking the supremum in n yields the second estimate in (A.3). 

For the proof of estimate (A. 4), take (f){x) = 1 on [n — + ^] and outside 



[n — 1, + 1], set w = (j)v, and compute similarly as the above. 

Using estimates in the above lemma, we can prove: 
Theorem A. 2 (local well-posedness in Hi). Suppose that 
(A.5) m"^ \\V\\l^ + \\V'\\l^ <oo. 

For any R> I, take 



□ 



T < mm{M-\R-^) 



we have 



(1) // ll^^oll//! < R, there exists a solution u{t) G C([0,T];if^) to (A.l) on [0,T] 
with initial data uq{x) satisfying 
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(2) This solution u{t) is unique among all solutions in C{[0,T]] H^) . 

(3) The data-to-solution map uq i— )■ u{t) is continuous as a mapping 
mT];Hl). 

Proof. We prove the existence by contraction in the space X, where 

X = {u \ \\u\\c{[0,T];H^) + \\Uxx\\i;jfL^Ll^^ H'^H'^n^'T 'f^Q„ - ^-^} ' 

where the constant C is chosen large enough to (10 times, say) exceed the imphcit 



constants in Lemma A.l Given u G X, let ^(u) denote the solution to 
(A.6) dMu) + dl^iu) - V^{u) = -2dx{u'') . 



with initial condition ip{u){0) = Uq. A fixed point ^{u) = -u in X will solve (A.l ). 



The local smoothing estimate (A. 3) applied to f = ^(u) and the estimate 



give the estimate 
(A.7) 



II </'(«) ||l2?l2 < Ikollffi +^l|M|li-m 5 



The maximal function estimate (A. 4) applied to f = ^(u) and the estimate 



{n%\\aLi<T'^'\ 



imply the estimate 

(A.8) \Mu)hlL^Ll ^\\uo\\L^,+Ty{u)\\L^Hi+T\\u\\l^Hi. 



Now applying to (A.6), and denoting v = (p{u)x instead: 

vt + Vxxx -Vv = -2{u^)xx + V'(^{u) . 



By (A. 3) again 



(A.9) y{u)x\\L^Li 



+ yiu)xx\\i^L2^Ll^ < holl^i + Uu^)xhlL2^Ll^ + ||I^VHIIl^L2 • 



Applying Gagliado-Nirenberg inequality to (f){x)u, where 0(x) = 1 on i'^ — ^, "^H- ^] and 
outside [n — l,n + 1], we obtain (writing Q for Qn and Q for Qn for the following): 



< i\\u\\r2 + \\u 



hence 



Win' 



ixWLi 



< 



\u. 



'x\\L 



2 \\U\ 
Q " ' 



< 



x\\L 



+ IkxIL?,) • 



Taking norm and applying the Holder inequality, we obtain 

\\iu^)x\\LlLl ^ \\Ux\\l^lI\\Al^L-^^{\\u\\lIlI + lkx||L2,L|) • 

Taking ^\ norm and applying the Holder inequality again yields 

\\{u^)xhiLlLl < l|Wa-||£-L-L2,l|M||^2L-L|(l|M||f2L2.L| + l^a- ||£2l2.L|) 
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Using the bounds || 



and 



we obtain 



|M||^2i2^2 < \\u\\L2 r2 < T^^lwllioo^a 



inserting into ( |A.9 ), 

(m)x||l5?L2 + \^{u)xx\t'^L%Ll, 



(A. 10) 



< 



Hi 



T ^Qn 



T OC T 2 



Summing (A.7), (A.8) and (A.IO), we obtain that ||v2(ti)IU < CR if \\u\\x < CR 



provided T < Comin(M ^,R with Cq small enough. Thus if : X — ?■ X. A 
similar argument establishes that <y9 is a contraction on X. 



Now suppose u,v e C{[0,T];H^) solve ( |A1] ). By (|X4|), 

(^•11) MM " / n n , n , ,|3 



vWilL^L^^ < Wvohl + T||t;||ioc^i + T\\v\ 



Set w = u — V. Then, with g = {u^ — v^) /{u — v) = u"^ + uv + f ^, we have 

Wt + W:rxx + 2(5'w)^ - Vw = . 



Apply (A. 3) to V = Wx, we obtain 

(A.12) Ww^Wl^lI < \\{9w).A\eiLlLl^ + \\V'w\\li,lI ■ 

The terms of ||(fi'w)2,'||^iL2L2 can be bounded in the following manner: 



(A. 13) 



UxVw\\iiL2^L2^ < \\Ux\\i^L^L^^^ \\vw\\iiL2 



< 



The term in the parentheses is bounded by 



which by (A. 11) and 
implies 



UxWe^L^L^^ < \\u\\l^hi 



\vhlLlLl^<T'^^\\v\\L^Ll 



UxVwy^LlLl^ ~ll«ll^oo^i,lbll,oc^i T^'WMiIl^lI+\Ml^hi)- 
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Same bounds follow for other terms in \\{gw)x\\ii l'^ l'^ , combined with ||V^'w||lil2 



< 



T\\V' \\l^\\w\\l^ Hii this estabhshes the estimate 
where the implicit constant depends on ||m||loo_h-i and ||f ||L°°_ffi- Same estimate follows 



for ||w||l^l2 by applying (A. 3) to v = w. Hence 



(A. 14) 



but applying (A. 4) to f = w yields 



(A. 15) 
since e.g. 



which can be proved by the same method as in (A. 13), and thus || (5''w)a;||j;^2 



< 



Substituting (A. 15) into (A. 14) implies w = for T sufficiently small, 



which then establishes the uniqueness of solutions in C([0,T]; H^). The continuity of 
the data-to-solution map can be proved by a similar argument. □ 

We now prove the global well-posedness in by (almost) conservation laws. 

Theorem A. 3 (global well-posedness). Suppose M < oo, where M is defined in 
(A. 5), for uq G , there is a unique global solution u G C^oc ( [0, C)o); if^) to (A.l) 



with \\U^L^Hi controlled by ||uo||f/i, T and M . 

Proof. First, note from Gaghado-Nirenberg inequality, ||'u||^4 ^ ||m|||2 ||Ma;||L2, we have 

\Wx\\\2 - \\u\\\2\\ux\\l^ < Hq{u) < ||Mx||i2 . 
Applying Peter-Paul inequality to the ||m|||2 H^ixlUa term gives us 



\u 



2 + \\ufj2 ~ Hoiu) + 



6 

L2 



(A.16) 



Suppose u solves (A.l), then 

jH,{u) = \{H',iu),JH',{u) + Vu)\ = \{H',iu),Vu)\ 

< M{\\Ux\\l2 + ||m||^2 + ||m||^4) < M{\\U.J12 + \\u\\l2 + ||M||i2||M:.||L2) 

< M(||«,,||i2 + ||n||i2 + \\u\\%) < M{H,{u) + ||n||i2 + Ml.) , 
on the other hand, by 

d 



dt 



P(u) = \{u,Vu)\<MP{u), 

and Gronwall inequality, we obtain a bound on in terms of ||mo||l2 and M, 

combine this with (A.16), and apply Gronwall again, we obtain a bound on Hq{u) 
and hence □ 
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Remark A. 4. A global well-posedness in for A; > 1 can in fact be proved, provided 
V & C^, by similar arguments. 
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